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Abstract: A measurement of the top-antitop (tt) charge asymmetry is presented using 
data corresponding to an integrated luminosity of 4.6 fb -1 of LHC pp collisions at a centre- 
of-mass energy of 7TeV collected by the ATLAS detector. Events with two charged leptons, 
at least two jets and large missing transverse momentum are selected. Two observables are 
studied: Aq based on the identified charged leptons, and Aq, based on the reconstructed 
tt final state. The asymmetries are measured to be 

Aq = 0.024 ± 0.015 (stat.) ± 0.009 (syst.), 

A% = 0.021 ± 0.025 (stat.) ± 0.017 (syst.). 

The measured values are in agreement with the Standard Model predictions. 
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1 Introduction 

The top quark is the heaviest elementary particle known to date. It was discovered in 1995 
at the Tevatron proton-antiproton (pp) collider by the CDF and DO collaborations [1, 2]. It 
is the only quark in the Standard Model (SM) that decays before hadronization occurs, and 
the only quark with Yukawa coupling to the Higgs boson close to unity. A precise study of 
top quark properties could shed light on possible physics models beyond the SM [3-9]. 

This analysis uses a data set corresponding to an integrated luminosity of 4.6 fb -1 of 
Large Hadron Collider (LHC) proton-proton (pp) collisions at a centre-of-mass energy of 
7TeV collected by the ATLAS detector. It is performed in the dilepton channel of the 
tt pair decay, realized when both W bosons decay to a charged lepton and a neutrino. 
The measured observables are the lepton-based charge asymmetry Aq and the tt charge 
asymmetry Aq. The observable Aq is defined as an asymmetry between positively and 
negatively charged leptons (electrons and muons) in the dilepton decays of the tt pairs, 

M _ N(A\ V \ > 0) - N(A\p\ < 0) 

C N(A\rj\ > 0) + A(A|r/| < 0) ’ 
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where 


A \v\ = \m+\ - \vi-\, (i-2) 

rj^+ ( 77 ^-) is the pseudorapidity 1 of the positively (negatively) charged lepton and N is 
the number of events with positive or negative A|ry|. While Aq is defined in eq. (1.1) 
as an asymmetry between positively and negatively charged lepton pseudorapidities, Aq 
corresponds to the asymmetry in top quark and antitop quark rapidities 2 , 

A tt = N(A\y\>0)-N(A\y\<0) 

C lV(A|y| > 0) + fV(A|y| < 0) ’ 

where 

A|y| = \yt\ ~ \Vi\, (1-4) 

yt {yi) is the rapidity of the top (antitop) quark, and N is the number of events with positive 
or negative A|y|. 

In SM tt production, the asymmetry is absent at leading-order (LO) in Quantum 
Chromodynamics (QCD) and is introduced by the next-to-leading-order (NLO) QCD con¬ 
tributions to the tt differential cross-sections, which are odd with respect to the exchange of 
t and t. At the LHC, the contributions to the asymmetries defined in eq. (1.1) and eq. (1.3) 
are predominantly from gg-initiated tt production, and gy-initiated production also has a 
non-negligible contribution. The yy-initiated processes are symmetric [10]. The asymmetry 
predicted in the SM is slightly positive, implying that the top quark is preferentially emitted 
in the direction of the quark in the initial state. In qq interactions at the LHC, the quark 
is in most cases a valence quark whereas the antiquark is from the sea. The asymmetry 
translates to a higher boost along the beam direction for the t-quark than for the f-quark. 
The rapidity distribution of the t is thus slightly broader than the one of the t. 

The SM predictions of Aq and Aq computed at NLO in QCD and including electroweak 
corrections (NLO QCD+EW) are [10] 


% = 0.0123 ± 0.0005 (scale), 

(1.5) 

# = 0.0070 ± 0.0003 (scale). 

( 1 . 6 ) 


These asymmetries are evaluated without acceptance cuts. The uncertainties are due to 
scale variations, estimated by simultaneous variation of the renormalization and factoriza¬ 
tion scale by a factor of half or two with respect to the reference scale value, which is 
set to the top quark mass. Recent next-to-next-to-leading order (NNLO) calculations of 
the forward-backward asymmetry for the Tevatron suggest that varying these scales signif¬ 
icantly underestimates the uncertainty due to higher order corrections [11], but no NNLO 
calculation has yet been published for pp interactions at the LHC energies. There is however 
a recent calculation obtained with the Principle of Maximum Conformality [12] which gives 
a consistent value of A^ = 0 . 01151 qqqo 3 (scale). The predicted value of Aq is smaller than 

1 The pseudorapidity is defined in terms of the polar angle 6 as y = — lntan(0/2). 

2 The rapidity is defined as y = \ In where E is the energy of the particle and p z is the component 

of the momentum along the LHC beam axis. 
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the prediction for Aq, since the directions of the leptons do not fully follow the direction 
of the parent t and t quarks. However, Aq can be measured more precisely, since it is 
determined without the need for a full reconstruction of t and t kinematics, which involves 
the use of jets and missing transverse momentum that are reconstructed with less precision 
than the kinematic variables of the leptons. The values of Aq and Aq as well as their 
correlation can be sensitive to new physics arising in top quark pair production [13-17]. 

The asymmetry Aq has been measured in the single-lepton decay channel by the AT¬ 
LAS [18] and CMS [19] collaborations at y/s = 7TeV. The CMS collaboration has reported 
measurements of Aq and Aq in the dilepton decay channel at yfs = 7TeV [20]. The mea¬ 
sured asymmetry values as well as those from a combination of ATLAS and CMS Aq results 
in the single-lepton decay channel [21] are consistent with the SM predictions. 

At the Tevatron collider, tt production has a forward-backward asymmetry with respect 
to the direction of the proton and antiproton beams. The asymmetry based on t and t 
quarks, A BB , is defined as 

A ti = N ( A U > 0) ~ N(Ay < 0) 

FB N(Ay > 0) + N(Ay < 0) ’ 

where 

Ay = y t -yt, (1-8) 

yt ( yi ) is the rapidity of the t (t) quark and N is the number of events with positive or 
negative Ay. An analogously defined lepton-based forward-backward asymmetry in tt pro¬ 
duction has been studied as well. At the Tevatron, tt events are predominantly produced by 
qq annihilation, thus the predicted asymmetries are typically larger than at the LHC, where 
gg-initiated production dominates. The Tevatron experiments have reported deviations of 
forward-backward asymmetries from the SM predictions [22, 23], which have motivated a 
number of further asymmetry measurements. Comparing the results with the latest NNLO 
calculations available at the Tevatron [11], the deviations reported by the CDF collabora¬ 
tion [24-26] are reduced, while the latest measurements by the DO collaboration [27, 28] 
are now in good agreement with the predictions. 

This paper is organized as follows. In section 2 the main components of the ATLAS 
detector relevant for this measurement are summarized. In section 3 the simulated samples 
used for the analysis are presented. In section 4 the object and event selection are described. 
In section 5 the kinematic reconstruction used for the Aq measurement is detailed. For 
comparison with theory prediction, the measurements are corrected for detector resolution 
and acceptance effects, as presented in section 6. Sections 7 and 8 describe the systematic 
uncertainties and the measurement results, respectively. Finally, the conclusions are given 
in section 9. 

2 The ATLAS detector 

The ATLAS detector [29] at the LHC covers nearly the entire solid angle around the colli¬ 
sion point." It consists of an inner tracking detector surrounded by a thin superconducting 

3 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in 
the centre of the detector and the 2 -axis along the beam pipe. The x-axis points from the IP to the centre 
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solenoid, electromagnetic and hadronic calorimeters, and a muon spectrometer incorporat¬ 
ing three large superconducting toroid magnets. The inner-detector system is immersed in 
a 2 T axial magnetic field and provides charged-particle-tracking in the range 1 77 I < 2.5. 

A high-granularity silicon pixel detector covers the interaction region and typically 
provides three measurements per track. It is surrounded by a silicon microstrip tracker 
designed to provide four two-dimensional measurement points per track. These silicon de¬ 
tectors are complemented by a transition radiation tracker, which enables radially extended 
track reconstruction up to 1 77 1 = 2.0. The transition radiation tracker also provides electron 
identification information based on the fraction of hits (typically 30 in total) exceeding an 
energy-deposit threshold corresponding to transition radiation. 

The calorimeter system covers the pseudorapidity range |? 7 | < 4.9. Within the region 
\r]\ < 3.2, electromagnetic calorimetry is provided by barrel and end-cap high-granularity 
lead/liquid-argon (LAr) electromagnetic calorimeters, with an additional thin LAr presam¬ 
pler covering \r]\ < 1.8 to correct for energy loss in the material upstream of the calorimeters. 
Hadronic calorimetry is provided by a steel/scintillator-tile calorimeter, segmented into 
three barrel structures within |? 7 | < 1.7, and two copper/LAr hadronic endcap calorime¬ 
ters. The solid angle coverage is completed with forward copper/LAr and tungsten/LAr 
calorimeters used for electromagnetic and hadronic measurements. 

The muon spectrometer comprises separate trigger and high-precision tracking cham¬ 
bers measuring the deflection of muons in a magnetic field generated by superconducting 
air-core toroids. The precision chamber system covers the region 1 77 I < 2.7 with three layers 
of monitored drift tube chambers, complemented by cathode strip chambers in the forward 
region. The muon trigger system covers the range |^| < 2.4 with resistive plate chambers 
in the barrel, and thin gap chambers in the endcap regions. 

A three-level trigger system is used to select interesting events. The Level-1 trigger is 
implemented in hardware and uses a subset of detector information to reduce the event rate 
to a design value of at most 75 kHz. This is followed by two software-based trigger levels, 
which together reduce the event rate to about 300 Hz. 

3 Simulated samples 

Several Monte Carlo (MC) simulated samples are used in the analysis to model the signal 
and background processes. The total background, estimated partly from these simulated 
samples, is subtracted from the data at a later stage of the analysis. The signal sample 
is used to correct the background subtracted data for detector, resolution and acceptance 
effects. The MC samples are also used to evaluate the systematic uncertainties of the 
measurement. 

The nominal simulated tt sample is generated using the PoWHEG-hvq [30-32] (patch4) 
generator with the CT10 [33] parton distribution function (PDF) set. The NLO QCD 
matrix element is used for the tt hard-scattering process. The parton showers (PS) and the 
underlying event are simulated using Pythia6 [34] (v6.425) with the CTEQ6L1 [35] PDF 

of the LHC ring, and the y-axis points upward. Cylindrical coordinates (r, (j>) are used in the transverse 
plane, (j> being the azimuthal angle around the beam pipe. 
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and the corresponding Perugia 2011C set of tunable parameters (tune) [36] intended to be 
used with this PDF. The hard-scattering process renormalization and factorization scales 
are fixed at the generator default value Q that is defined by 

Q = +p\, (3.1) 

where mt and pt are the top quark mass and the top quark transverse momentum, evaluated 
for the underlying Born configuration (i.e. before radiation). Additional tt samples used 
to evaluate signal modelling uncertainties are described in section 7. Signal samples are 
normalized to a reference value of cr t j = 177^^ pb for a top quark mass of mt = 172.5 GeV. 
The cross-section has been calculated at next-to-next-to-leading-order (NNLO) in QCD 
including resunmiation of next-to-next-to-leading logarithmic (NNLL) soft gluon terms [37— 
42] with top++2.0 [43]. The PDF and strong coupling (ct s ) uncertainties were calculated 
using the PDF4LHC prescription [44] with the MSTW2008 68 % CL NNLO [45, 46], CT10 
NNLO [33, 47] and NNPDF2.3 5f FFN [48] PDF sets, and added in quadrature to the 
scale uncertainty. The NNLO+NNLL cross-section value is about 3% larger than the exact 
NNLO prediction, as implemented in Hathor 1.5 [49]. 

The MC generators which are utilized to estimate the backgrounds are as follows. 
Single-top processes in the Wt channel are generated with the MC@NLO event generator 
(v4.01) [50, 51] with the CT10 PDF. The parton showers, hadronization and the underlying 
event are modelled using the HERWIG (v6.520) [52, 53] and JIMMY (v4.31) [54] generators. 
The CT10 PDF with the corresponding ATLAS AUET2 tune [55] is used for parton shower 
and hadronization settings. For Z/ 7 *-fjets and diboson events ( WW , WZ and ZZ ), Alp- 
GEN (v2.13) [56] interfaced to Herwig and Jimmy is used. The CTEQ6L1 PDF and the cor¬ 
responding ATLAS AUET2 tune is used for the matrix element and parton shower settings. 
The Wt background process is normalized to the reference NLO+NNLL QCD [57] predic¬ 
tion. Diboson production is normalized to the reference NLO QCD prediction obtained 
using MCFM [58] and MC@NLO generators with the MSTW2008 NLO PDF [45]. The 
Z/ 7 * — > ee//z/i+jets cross-section is normalized using a control region in data as detailed 
in section 4. The Z/ 7 * —> rr+jets events are normalized to a NNLO reference cross-section 
using the FEWZ [59] and ZWPROD [60] programs with the MSTW2008 NNLO PDF. 

To realistically model the data, the simulated samples are generated with an average 
of eight additional inelastic pp interactions from the same bunch crossing (referred to as 
pileup) overlaid on the hard-scatter event. Simulated samples are processed through AT¬ 
LAS detector simulation. For the majority of the samples, a full detector simulation [61] 
based on GEANT4 [62] is used. Some of the samples used for assessment of generator 
modelling uncertainties are obtained using a faster detector simulation program that relies 
on parameterized showers in the calorimeters [61, 63]. Simulated events are then processed 
using the same reconstruction algorithms and analysis chain as the data. 

4 Object and event selection 

The data sample collected by the ATLAS detector in 2011 at a centre-of-mass energy of 
7TeV is used for the analysis. The integrated luminosity of the sample is 4.6 fb -1 with an 
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overall uncertainty of 1.8% [64], The analysis makes use of reconstructed electrons, muons, 
jets and missing transverse momentum in the detector. Electrons are reconstructed as clus¬ 
ters of energy deposits in the electromagnetic calorimeter, matched to a track in the inner 
detector. They are required to pass a set of tight selection criteria [65]. The selected elec¬ 
trons have to satisfy a requirement on their transverse energy (E?x) and the pseudorapidity 
of the associated calorimeter cluster (|r? c iuster|) : -Et > 25GeV and |r? c i us ter| < 2.47. The 
electrons in the region 1.37 < |7? c i U ster| < 1-52, which corresponds to a transition between 
the barrel and endcap electromagnetic calorimeters, are excluded. Electrons are required to 
be isolated, using the requirements described as follows (excluding calorimeter deposits and 
tracks from the electrons). The E x within a cone of size A R = yj (Ar/) 2 + (A cj)) 2 = 0.2 and 
the scalar sum of track px within a cone of A R = 0.3 around the electron are required to 
be below E^- and r/-dependent thresholds. The efficiency of this isolation requirement on 
electrons is 90%, and its goal is to reduce the contribution from hadrons mimicking lepton 
signatures, as well as leptons produced in heavy-hadron decays or photon conversion. These 
are referred to as fake and non-prompt leptons (NP) in the following. 

Muons are reconstructed by matching a track in the inner detector to a track segment 
in the muon spectrometer. They are required to pass tight selections [66]. The selected 
muons are required to have px > 20GeV and 1 77 I < 2.5. To reject fake and non-prompt 
muons, the following isolation requirements are imposed: the calorimeter transverse energy 
within a cone of A R = 0.2 around the muon is required to be less than 4GeV and the scalar 
sum of track px within a cone of A R = 0.3 is required to be less than 2.5 GeV (excluding 
the calorimeter deposits and tracks from the muons). 

Jets are reconstructed from energy deposits in the calorimeter, using the anti-Ap al¬ 
gorithm with a distance parameter R = 0.4 [67]. The energy of the input clusters [68] is 
corrected to the level of stable particles using calibration factors derived from simulation 
and data [69]. The jets are required to have a px of at least 25 GeV and |p| < 2.5. To 
suppress the contribution from low-px jets originating from pileup interactions, tracks as¬ 
sociated with the jet and emerging from the primary vertex are required to account for at 
least 75% of the scalar sum of the px of all tracks associated with the jet. A primary ver¬ 
tex, originating from pp interactions, is a reconstructed vertex required to have at least five 
associated tracks with px > 0.4 GeV. In the cases where more than one primary vertex is 
reconstructed, the vertex with the highest Xark^T chosen and assumed to be associated 
with the hard-process, and the sum runs over all associated tracks. 

The missing transverse momentum (_E™ 1SS ) is a measure of transverse momentum im¬ 
balance due to the presence of neutrinos. It is reconstructed from the transverse momenta of 
jets in the kinematic range of px > 20 GeV and \rj\ < 4.5, electrons, muons, and calorimeter 
clusters not associated with any of the reconstructed objects, as detailed in ref. [70]. 

Using the objects reconstructed as above, an event selection optimized for signatures 
corresponding to tt events in which both W bosons from the t and t quarks decay to leptons 
is performed. Events are required to have been selected by a single-electron trigger with a 
threshold of 20 or 22GeV (depending on the data-taking period), or a single-muon trigger 
with a threshold of 18 GeV. They are required to have exactly two isolated, oppositely 
charged, leptons. Depending on the lepton flavours, the sample is divided into three analysis 
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channels referred to as ee, e/x and /x/x. To reduce the Drell-Yan production of Z/ 7 *+jets 
background, the invariant mass of the two leptons (mu) is required to be above a threshold 
used to suppress 7 * —> it production background and outside a Z boson mass window in 
the ee and /x/x channel events. The following requirements are used: mu > 15 GeV and 
\mu~mz\ > 10 GeV. In the ee and /x/x channels the Drell-Yan and diboson backgrounds are 
further reduced using a requirement on the missing transverse momentum, E!f lss > 60 GeV. 
In the e/x channel the Z/ 7 *+jets background is smaller and suppressed by requiring the 
scalar sum of the p t of the two leading jets and leptons (FTp) to be larger than 130 GeV. 

The background contributions are estimated using a combination of techniques us¬ 
ing data and Monte Carlo events. In the case of single-top and diboson processes, both 
the shape and normalization of the distributions are taken from the simulation. For 
Z/ 7 * —> ee//x/x+jets events, simulated MC events are used to model the shape of the 
distributions, but a data control region is used for normalization. Drell-Yan events with 
^nnss > go GeV are affected by energy mismeasurements, that are difficult to model in 
simulation. A control region with events with mu in the Z-rnass region is defined to study 
the effect of mismeasured E!f lss . The relative E™ 1SS , defined as the projection of the missing 
transverse momentum onto the direction of the jet or charged lepton with closest <f), is used 
to identify the events with mismeasured objects. Events with energy mismeasurements are 
characterized by high values of relative E™ ss . A cut is applied to the relative _E™ 1SS , and 
data and simulation are then compared to derive a normalization correction factor which is 
applied to the simulated sample. The Z/ 7 * —>• rr contribution is estimated from MC simu¬ 
lation. The background stemming from events with at least one non-prompt or fake lepton 
is estimated from the data, since the lepton misidentihcation rates are difficult to model 
in MC simulation. A matrix method technique is used [71]. It consists of selecting data 
samples dominated either by real leptons or by fake leptons, and estimating the efficiencies 
for a real or fake lepton to satisfy the isolation criteria. 

After the final selection, the data sample contains more than 8000 events, with an 
expected signal-to-background ratio of approximately six. The number of events in data and 
simulation, including statistical and systematic uncertainties, are compared in table 1. After 
selection, the largest number of events is observed in the e/x channel, which has the highest 
branching ratio and the loosest background suppression cuts. The ee channel has the lowest 
number of events because of the stringent requirements on lepton kinematics. Figure 1 shows 
good agreement between the data and the SM predictions for the jet multiplicity, lepton 
Pt an d lepton pseudorapidity distributions. The A|ry| distributions are shown in figure 2 
for the three channels separately. 

5 Kinematic reconstruction 

For the measurement of the tt charge asymmetry, the direction of the top and antitop quarks 
needs to be determined. The four-momenta of top quarks in selected events are computed 
with a kinematic reconstruction using the objects observed in the detector. The reconstruc¬ 
tion is based on solving the kinematic equations obtained when imposing energy-momentum 
conservation at each of the decay vertices of the process. In the dilepton channel, at least 
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Channel 

ee 

efi 


tt 

621 

± 

5 

± 

59 

4670 

± 

10 

± 

325 

1780 

± 

10 

± 

120 

Single top 

31.6 

± 

1.7 

± 

3.8 

230 

± 

5 

± 

21 

83.9 

± 

2.7 

± 

8.3 

Diboson 

22.8 

± 

0.9 

± 

2.6 

177 

± 

3 

± 

16 

61.5 

± 

1.5 

± 

6.1 

Z — * ee (DD) 

20.8 

± 

1.7 

± 

1.4 



— 





— 



Z —> fifi (DD) 



— 



2.1 

± 

0.5 

± 

0.7 

77 

± 

4 

± 

12 

Z — y tt 

18.6 

± 

1.8 

± 

7.0 

170 

± 

6 

± 

60 

67 

± 

4 

± 

25 

NP & fake (DD) 

19 

± 

4 

± 

19 

99 

± 

10 

± 

63 

26.8 

± 

5.1 

± 

1.9 

Total expected 

734 

± 8 

± 63 

5350 

± 20 

±340 

2100 

± 10 

±130 

Data 

740 

5328 

2057 


Table 1: Observed number of data events in comparison to the expected number of signal 
events and all relevant background contributions after the event selection. The backgrounds 
are estimated from the MC simulation or from the data-driven methods (DD) described in 
section 4. Events with one or more non-prompt or fake leptons are referred to as "NP &; 
fake". The first uncertainty is statistical, the second corresponds to systematic uncertainties 
on background normalization and detector modelling described in section 7. The values 
labeled with "—" are estimated to be smaller than 0.5. 


two neutrinos are produced and escape undetected. Consequently, the system is undercon¬ 
strained and its kinematics cannot be fully determined without further assumptions (for 
example on the W boson and top quark masses, and the pseudorapidities of the neutrinos 
from the W boson decays). Moreover, several ambiguities have to be resolved to find the 
correct solution. For example, the lepton and jet from the same decay chain have to be 
associated. In an event with two leptons and two jets, this leads to two possible associa¬ 
tions. In this analysis, the neutrino weighting technique [72] is used. This procedure steps 
through different hypotheses for the pseudorapidity of the two neutrinos in the final state. 
These hypotheses are made independently for the two neutrinos. For each hypothesis, the 
algorithm calculates the full event kinematics, assuming the W boson and the top quark 
masses. It then assigns a weight to the resulting solution based on the level of agreement 
between the calculated and measured missing transverse momentum. The weight is defined 
as 


w = 


II exp 
d=x,y 




miss,calc 


-E 


2 ( 


(?i 


miss,obs 
d _ 

2 



(5.1) 


with £ , “ uss ’ ob& being the projection of the measured missing transverse momentum along the 
axes defining the transverse plane (d = x, y) and £^ mss ’ calc the projection calculated with the 
assumed r/ values of the neutrino pair. The resolution on the missing transverse momentum 
is denoted a E mi SS , and defined as (7 E mi SS = 0.5y^2Er GeV [70]. The total transverse energy, 
^ -&T is defined as E t = E{ sin 9i where E{ and 9{ are the energy and the polar 

angle of calorimeter cells associated with clusters. All possible lepton-jet associations are 
considered and jet energy mismeasurements are accounted for by random shifts of the jet 
energies within their resolutions. The solution corresponding to the maximum weight is 
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Figure 1: Comparison of the expected and observed distributions of (a) the jet multiplicity, 
(b) the lepton transverse momentum p^, (c) the lepton pseudorapidity r] and (d) the missing 
transverse momentum F™ 1SS , shown for the combined ee, ep and /jfi channels. Events 
beyond the range of the horizontal axis of (a), (b) and (d) are included in the last bin. The 
hatched area corresponds to the combined statistical and systematic uncertainties. Events 
with one or more non-prompt or fake leptons are referred to as "NP & fake". 


selected to represent the event. 

As a result of the scan over neutrino pseudorapidities and jet energy values, the recon- 
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Figure 2: Comparison of the expected and observed distributions of the A\r)\ variable for 
the (a) ee, (b) efj, and (c) /r/r channels. The hatched area corresponds to the combined sta¬ 
tistical and systematic uncertainties. Events with one or more non-prompt or fake leptons 
are referred to as "NP & fake". 


struction efficiency, corresponding to the fraction of events in which solutions for t and t 
four-momenta are found, is estimated to be about 80% in the data. In the other 20% of 
events, no solution to the system of the kinematic equations could be found, and the events 
are not used for the measurement of Aq. The performance of the reconstruction algorithm 
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for key variables, such as the top quark rapidities and A|y|, is evaluated using the nominal 
tt simulated sample. The fraction of reconstructed MC events where the sign of A|y| is 
determined correctly is about 70%. In the simulated samples, the correct combination of 
the charged leptons and two jets from 6(6)-quarks is found in approximately 80% of the 
events with exactly two reconstructed jets, both of which are matched to the 6(6)-quarks. In 
case all events passing the event selection are considered, the correct combination is found 
in approximately 47% of the events. 

In figure 3 the distributions of the top quark transverse momentum, top quark rapidity 
and the tt invariant mass are shown for the combined ee, e/i and /ijU channels. In figure 4 
the A|y| distribution is shown separately for each of the ee, e/r and /i/r channels. Good 
agreement between the observed and expected distributions is found. 
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Figure 3: Comparison of the expected and observed distributions of (a) the top and 
antitop quark transverse momentum p x, (b) top and antitop quark rapidity and (c) the 
tt invariant mass, shown for the combined ee, ep and pp channels. The hatched area 
corresponds to the combined statistical and systematic uncertainties. Events with one or 
more non-prompt or fake leptons are referred to as 11 NP &; fake". 
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Figure 4: Comparison of the expected and observed distributions of the A|y| variable for 
the (a) ee, (b) e^i and (c) /r/i channels. The hatched area corresponds to the combined sta¬ 
tistical and systematic uncertainties. Events with one or more non-prompt or fake leptons 
are referred to as "NP &; fake". 


- 13 
















































ATLAS Simulation, e|x, is = 7 TeV 


1 

- 0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

3 

94 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

92 

6 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

93 

5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

93 

4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

93 

4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

93 

4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

94 

4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 

94 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 

93 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 

93 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 

93 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

6 

93 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

5 

92 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

- 95 

_1_ 

2 

0 

0 

0 

0 

0 

_ 

0 

0 

0 

0 

0 

0 

_ 

0 


■_i_i_i_i_i_i_i_i_i_ 

-2 0 2 

Generated A|r|| 


ATLAS Simulation, e|r, is = 7 TeV 


1 i T T 1 1 1 ^ 

1 1 1 

1 1 l 1 1 1 l 1 

; 3 

6 16 

57 

12 

28 50 

29 

29 

49 28 

11 

56 

, 1 1 1 1 1 , 1 

17 6 

1 1 . 

3 

1 1 1 1 1 1 1 1 

-4 -2 

0 

2 4 


Generated A|y| 


(b) 


Figure 5: Response matrices for (a) the lepton A|? 7 | and (b) tt A\y\ observables in the 
e/i channel. Each column of the matrices is normalized to unity and values are reported as 
percentage (%) units. Values smaller than 0.5% are rounded to 0%. 


6 Corrections 

For comparison with theoretical calculations, the measurements are corrected for detector 
resolution and acceptance effects. The corrections are applied to the observed A|??| and A|y| 
spectra. Apart from the corrected inclusive asymmetry values, particle- or parton-level A\r]\ 
and A|y| distributions are obtained and presented as normalized differential cross-sections 
in section 8 . Acceptance corrections are included, thus all the results correspond to an 
extrapolation to the full phase-space for tt production. 

In case of the the resolution of the measured lepton A|r/| is very good. Figure 5(a) 
shows, for the ey channel, the probability of an event with a generated value A|^| in the j-th 
bin to be reconstructed in the z-th bin of the corresponding distribution. This probability 
distribution is defined to be the response matrix for the observable A1 77 1. The diagonal bins 
of the response matrix account for more than 90% of the events. The acceptance and the 
small migrations are accounted for by the bin-by-bin correction described in subsection 6 . 1 . 

In case of the A the top quark direction, which is necessary to determine the tt 
asymmetry, is evaluated using the kinematic reconstruction of the events, described in 
section 5. In addition to lepton directions and energies measured with very good resolution, 
jet four-momenta and F™ 1SS measured with worse resolution are used in reconstructing the 
t and t four-momenta. The resolution of tt A|y| (figure 5(b)) is thus much worse than that 
for the lepton A|? 7 |. In order to correct for detector resolution and acceptance effects in 
the fully Bayesian unfolding (FBU) technique [73] described in subsection 6.2 is used. 

6.1 Correction of the lepton-based asymmetry 

For Aq , bin-by-bin correction factors that also extrapolate to the full acceptance for the 
tt production are used. The goal of this procedure is to find an estimate of the true 
distribution, given an observed distribution and an expected background distribution. For 
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the lepton-based results, true distributions are obtained at particle level using leptons before 
Quantum Electrodynamics (QED) final-state radiation 1 . The following notation is used: /z 
and p are vectors of true distribution values and its estimate, respectively. An observed 
distribution is denoted by n and its expected value from simulation by i/ MC . An expected 
background distribution is denoted by (3. For the z-th bin of the asymmetry distribution, 
the estimate of the true value is obtained by applying a correction factor Ci to the difference 
between the observed number of events and the expected number of background events, 


pi — Ciijli f3i) . 


( 6 . 1 ) 


The Ci are estimated using the tt MC simulated sample as 


Ci = 


7 /MC : 


( 6 . 2 ) 


where and uf lc are the predictions for the number of events in the z-th bin of the true 
and reconstruction-level distributions. 

The bin-by-bin correction of Aq is tested on simulation samples reweighted such that 
different levels of asymmetry A|zy| are introduced. Samples are reweighted according to a 
linear function of A\r]\ with a slope between —6% and 6% in steps of 2%. Corrected values 
obtained from reweighted distributions are found to be in good agreement with the input 
value, following a linear relationship. The choice of the binning is done by optimizing the 
linearity of the method and the expected statistical uncertainty of the asymmetry. The 
results in section 8 are obtained with A|?y| distribution binned in 14 bins in the interval 
between —3 and 3. 

The correction factors depend strongly on the channel and the bin, with the outer bins 
receiving larger fractional corrections. The ee channel has the lowest acceptance and thus 
the highest correction factors, reaching values of 500 in the outer bins, in which the events 
are mostly outside the detector fiducial acceptance. The e/z channel has a much higher 
acceptance, and the correction factors vary between 10 and 60. The dependence of the 
correction factors on the MC model and PDF is small, up to approximately 5%. 


6.2 Unfolding of the tt asymmetry 

In case of sizeable migrations across the bins of the considered distribution, the migrations 
need to be taken into account without introducing a significant bias during the correction 
procedure. Unfolding is better suited for the purpose than the bin-by-bin correction factors 
described in subsection 6.1. Using the response matrix ( R ), the true distribution (/z) is 
related to the expected reconstruction-level distribution (iz) and the expected background 

(0) by 

v = R/z + (3. (6-3) 

4 The particle-level definition uses status-code 3 for Pythia6 for electrons and muons produced in W 
boson decays. In addition, electrons and muons produced from status-code 3 r leptons are used. These 
particles are used both for the unfolding and for the predictions of the MC generators. 
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In the FBU technique, the maximum likelihood estimator of p. L(/h), is given by 

N 

log L(n) = ^log P{ni\Vi) -aS{p) ; p(p) oc L(p) , (6.4) 

2—1 

with P the Poisson distribution, n the observed distribution, S a regularization function and 
a a regularization parameter. The sum in i runs over all N bins of the distributions. The 
probability density of the unfolded spectra p(p) is proportional to L(p). The regularization 
function S is selected such that the spectra with a desired quality, such as smoothness, are 
preferred. The regularization parameter a controls the relative strength of the regularization 
when evaluating the likelihood. The unfolded spectrum and its associated uncertainty are 
extracted from the probability density p(p). The statistical uncertainty corresponds to 
the width of the shortest interval covering 68 % probability, and the unfolded spectrum 
corresponds to the middle of that interval. 

The response matrix is obtained using information from the nominal tt simulated sam¬ 
ple and, in particular, using the top quarks before their decay (parton level) and after QCD 
radiation 5 . 

As explained for the lepton-based asymmetry, the correction is done at the level of true 
dilepton events (where the two top quarks decay to electrons or muons, either from a direct 
W boson decay or through an intermediate r lepton decay). 

Using the vector of the true distribution’s estimated values p, the regularization func¬ 
tion is defined based on the curvature S(iu) = \C(p) — C(p)\, with 

JV—l 

C\p) — ' (Aj+i,i , and Aj_|_pj — Pi+i Ah- (6.5) 

i=2 

As in the case of the lepton-based asymmetry, linearity tests are performed. A given 
asymmetry value is introduced by reweighting the samples according to a linear function of 
tt A|y| with a slope between - 6 % and 6 % in steps of 2%. Unfolded values obtained from 
reweighted distributions are observed to be in good agreement with the injected values 
of Aq, following a linear relationship. This linearity test is performed with and without 
regularization and yields similar performance. The binning used for the A|y| distribution as 
well as the regularization parameter are optimized simultaneously to minimise the expected 
statistical uncertainty while achieving good linearity. The results in section 8 are obtained 
with a regularization parameter a = 10”'. The A|y| distribution is binned in 4 bins in the 
interval between —5 and 5. For this binning choice, at least 50% of the events populate the 
response matrix diagonal bins for each of the ee, ep and pp channels (figure 5(b)). 

The overall correction which is applied to the distribution varies between factors of 10 
and 100, depending on the channel and the bin. As shown in figure 5 the bins used for 
the tt A|y| distribution are wider than the bins used for the lepton A1 77 1 distribution. The 
acceptance correction applied to the outer bins of A|y| is thus smaller than the correction 
obtained for the outer bins of A|r;| distribution. 

5 The parton-level definition uses status-code 155 for Herwig and 3 for Pythia6 for both the unfolding 
and for the predictions of the MC generators. 
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ee 

e/r 

PP 

comb. 

Measured value 
Statistical uncertainty 

0.101 

±0.052 

0.009 

±0.019 

0.047 

±0.030 

0.024 

±0.015 

Lepton reconstruction 

±0.011 

±0.008 

±0.009 

±0.008 

Jet reconstruction 

±0.006 

±0.001 

±0.004 

±0.001 

ZT'miss 

TZ/rj-i 

±0.001 

<0.001 

±0.002 

<0.001 

Signal modelling 

±0.004 

±0.003 

±0.003 

±0.003 

PDF 

±0.004 

<0.001 

<0.001 

<0.001 

NP &; fake 

±0.016 

±0.001 

±0.001 

±0.001 

Background 

±0.003 

±0.002 

<0.001 

±0.001 

Total sys. 

±0.021 

±0.009 

±0.012 

±0.009 


Table 2: Measured value and uncertainties for the lepton-based asymmetry Aq. Uncer¬ 
tainties with absolute value below 0.001 are considered negligible for the total uncertainty. 

7 Systematic uncertainties 

The uncertainties of the Aq corrections and A ^ unfolding method are estimated from the 
non-closure in the linearity test described in section 6. For Aq the uncertainties are —0.002 
in ee channel and negligible (< 0.001) in the e/r and /r/i channels. For Aq the uncertainties 
are 0.002 in the ee channel and negligible (< 0.001) in the e^u and pp, channels. For both 
Aq and Aq the uncertainties have a negligible contribution to the measurement uncertainty 
and are not considered for the evaluation of the total systematic uncertainty of the results. 

The systematic uncertainties considered in this analysis are classified into three cate¬ 
gories: detector modelling uncertainties, signal modelling uncertainties and uncertainties 
related to the estimation of the backgrounds. The contributions of these sources of un¬ 
certainty are summarized in table 2 for the lepton-based asymmetry Aq and in table 3 
for the tt asymmetry A The resulting variations are assumed to be of the same size in 
both directions and are therefore symmetrized. Apart from one-sided uncertainties, as in 
the case of the comparison of different MC models, the symmetrization does not notably 
modify the uncertainty values. 

Detector modelling uncertainties are evaluated by performing corrections for detec¬ 
tor effects for Aq and A^, with the response matrices corresponding to the systematic 
variations. Effects of detector modelling uncertainties on the background are included by 
subtracting the background, varied accordingly, from the data. The following sources are 
considered. 

• Lepton reconstruction 

The uncertainty due to lepton reconstruction includes several sources. Lepton mo¬ 
mentum scale and resolution modelling correction factors and associated uncertainties 
are derived from comparisons of data and simulation in Z -^r ti events [65, 66]. Uncer- 
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ee 

e/i 

IJA 

comb. 

Measured value 
Statistical uncertainty 

0.025 

±0.069 

0.007 

±0.032 

0.043 

±0.045 

0.021 

±0.025 

Lepton reconstruction 

±0.008 

±0.008 

±0.004 

±0.007 

Jet reconstruction 

±0.015 

±0.009 

±0.006 

±0.009 

ZT'miss 

TZ/rj-i 

±0.015 

±0.005 

±0.008 

±0.007 

Signal modelling 

±0.004 

±0.004 

±0.003 

±0.003 

PDF 

±0.004 

±0.005 

±0.004 

±0.005 

NP Sz fake 

±0.013 

±0.011 

±0.003 

±0.008 

Background 

±0.003 

±0.002 

±0.004 

±0.003 

Total sys. 

±0.027 

±0.018 

±0.013 

±0.017 


Table 3: Measured value and uncertainties for the tt asymmetry Aq. 

tainties in the modelling of trigger, reconstruction and lepton identification efficiencies 
are also included. Data-to-simulation efficiency corrections, and their uncertainties, 
are derived from J/ij) —>• tt, Z —> tt and W —> ev events. 

• Jet reconstruction 

The effects include the jet energy scale and jet resolution uncertainties. Jet energy 
scale uncertainty is derived using information from test-beam data, LHC collision 
data and simulation [69]. It includes uncertainties in the flavour composition of the 
samples and mismeasurements due to close-by jets and pileup effects. Jet energy 
resolution and reconstruction efficiency uncertainties are obtained using minimum 
bias and QCD dijet events [69, 74]. 

• E™ iss 

The uncertainties from the energy scale and resolution corrections for leptons and 
jets are propagated to the E™ 1SS . The category accounts for uncertainties in the 
energies of calorimeter cells not associated with the reconstructed objects and the 
uncertainties from cells associated with low-px jets (7 GeV< pt < 20 GeV) [70] as 
well as the dependence of their energy on the number of pileup interactions. 

The uncertainties due to the modelling of the signal tt distributions are evaluated by 
performing the linearity test for signal model samples generated with various assumptions. 
The following sources are quoted. 

• Signal modelling 

The uncertainty is evaluated by adding in quadrature the MC generator uncertainties, 
initial- and final-state radiation (ISR and FSR), underlying event (UE) and colour re¬ 
connection (CR) uncertainties described in the following. The systematic uncertainty 
related to the choice of a MC generator includes the difference between the nom¬ 
inal sample generated with PoWHEG-hvq + PYTHIA6 and samples generated with 
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MC@NLO + Herwig, PowHEG-hvq + Herwig and Alpgen + Herwig. The ef¬ 
fects of renormalization and factorization scale choice are evaluated with a dedicated 
pair of samples generated with MC@NLO + HERWIG. In these samples renormal¬ 
ization and factorization scales are varied simultaneously by a factor of two with 
respect to the reference scale. The reference scale is fixed at the MC@NLO gener¬ 
ator default, which is defined as the average of the t and the t transverse masses, 
Q = y // l/2(pxt +PTf) + where PTt(t) corresponds to the transverse momentum 
of the t or t. Since the effects covered by generator comparisons and scale variations 
partially overlap, only the largest contribution from all comparisons is used. For the 
lepton-based asymmetry the dominant contribution was found to stem from the differ¬ 
ence between the nominal sample and the sample generated with ALPGEN + HERWIG. 
For the ti asymmetry the contributions of the comparison of the baseline sample re¬ 
sult to the results obtained with each of MC@NLO + Herwig, PowHEG-hvq + 
Herwig and Alpgen + Herwig samples are of comparable size and significantly 
larger than the contribution of the renormalization and factorization scale uncertainty. 
The amount of ISR and FSR are treated as an additional source of signal modelling 
uncertainty. It is evaluated using samples generated with ALPGEN + PYTHIA6 with 
variations of parameters controlling the renormalization scale used in ALPGEN and in 
the PYTHIA6 parton shower. The renormalization scale is varied by factors of 0.5 and 
2. The PYTHIA6 settings correspond to Perugia radLO and radHi tunes [36]. Apart 
from this, the UE and CR uncertainties are evaluated by comparing samples gen¬ 
erated with PoWHEG-hvq + PYTHIA6, using Perugia2011, Perugia2011 mpiHi and 
Perugia2011 noCR tunes [36]. For Aq . the contributions from the choice of MC gen¬ 
erator and from ISR and FSR exceed the non-perturbative UE and CR contributions. 
For Aq, the contributions from the choice of MC generator and from radiation and 
non-perturbative modelling uncertainties are comparable. 

• PDF uncertainty 

The uncertainty due to the PDF is evaluated by performing linearity tests with sam¬ 
ples obtained from the nominal signal sample, generated with CT10 PDF, reweighted 
to other PDFs. The CT10 error set as well as MSTW2008 68% CL NLO [45] and 
NNPDF2.3 NLO (ct s = 0.118) [48] central predictions are used. For each asymmetry 
value, the largest value of the three sources is quoted as uncertainty. 

The uncertainties on the modelling of the SM backgrounds are divided into two cate¬ 
gories described below. 

• NP & fake 

This source corresponds to the uncertainty in the estimation of processes fulfilling 
the event selection due to non-prompt or misidentihed leptons. The uncertainties are 
obtained by varying the efficiencies for a real or fake lepton to pass the tight selection, 
and are affecting both the normalization of the background and its shape. 

• Background 

The uncertainties in the modelling of diboson, Z+jets and single-top SM processes are 
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quoted in the background category. They are evaluated by varying the normalization 
of each of these processes by the uncertainty on its cross-section. The uncertainty on 
the overall luminosity of 1.8% is also entering this category [64], 

For both the lepton-based asymmetry Aq and the tt asymmetry Aq, the statistical 
uncertainty is larger than the total systematic uncertainty. The Aq measurement has a 
combined statistical uncertainty of 1.5%, whereas the combined systematic uncertainty is 
0.9%. The largest source of Aq systematic uncertainty is the lepton reconstruction un¬ 
certainty, which accounts for approximately 90% of the total systematic uncertainty. The 
uncertainty on the asymmetry Aq measured in the ee channel receives a sizeable contri¬ 
bution from the NP & fake leptons category (1.6%). This, however, does not significantly 
impact the combined systematic uncertainty since the ee channel receives a small weight in 
the combination, as detailed in section 8 . The tt asymmetry Aq has a combined statistical 
uncertainty of 2.5% and a combined systematic uncertainty of 1.7%. The detector modelling 
uncertainties account for approximately 80% of the combined systematic uncertainty, with 
comparable large contributions from the lepton reconstruction, the E™ 1SS (0.7%) and the 
jet reconstruction uncertainty (0.9%). The NP & fake contribution to the Aq systematic 
uncertainty is also sizeable ( 0 . 8 %). 

The uncertainties related to detector and background modelling are evaluated in each 
bin of the corrected distributions and presented in section 8 . 

8 Results 

After the event selection and reconstruction but before the correction described in section 6 
the inclusive lepton and tt asymmetries measured in the data are Aq — 0.021 ± 0.011 
(stat.) and Aq — 0.003 ± 0.012 (stat.), respectively for the combination of the ee,e/j, and 
/ ifj, channels. After the subtraction of the background contribution, the measured data 
asymmetries are Aq — 0.029 ± 0.013 (stat.) and Aq = 0.006 ± 0.014 (stat.). The corre¬ 
sponding asymmetry predictions in the nominal simulated tt sample are Aq — 0.005 ± 0.003 
(stat.) and Aq = 0.008 ± 0.003 (stat.). This sample is generated with the PoWHEG-hvq + 
PYTHIA6 generator with a particle-level lepton asymmetry of Aq — 0.0045 ±0.0009 (stat.) 
and a parton-level tt asymmetry of Aq — 0.0071 ± 0.0009 (stat.), evaluated in the full 
phase-space. 

After the correction for detector, resolution and acceptance effects, the normalized 
differential cross-sections corrected to particle and parton level are obtained for A|iy| and 
A|y| separately for the three channels. From these distributions, the inclusive asymmetry 
values can be extracted. The inclusive results obtained in the ee, e/i and /i/U channels 
(see tables 2 and 3) are then combined using the best linear unbiased estimator (BLUE) 
method [75, 76]. All systematic uncertainties are assumed to be 100% correlated, except 
for the uncertainties on electrons and muons and on the NP &; fake lepton background. 

The normalized differential cross-sections for A| 77 I and A|y| are presented in figure 6 
for the efj, channel. Good agreement is observed between the measured distributions and 
the ones predicted by PoWHEG-hvq + PYTHIA6. The normalized differential cross-sections 
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Figure 6: Normalized differential cross-sections for (a) lepton A\r]\ and (b) tt A\y\ in the 
ep channel after correcting for detector effects. The distributions predicted by Powheg- 
hvq + PYTHIA6 are compared to the data in the top panels. The bottom panels show the 
ratio of the corrected data to the predictions. The error bars correspond to the statistical 
uncertainties and the hatched area to the systematic uncertainties. 


in that channel are also presented with statistical and systematic uncertainties in tables 4 
and 5. The systematic uncertainties for the differential distributions do not include the 
signal modelling uncertainties, which could not be evaluated with sufficient precision due 
to the limited statistics of the simulated samples. For both distributions, the statistical 
uncertainty is somewhat larger than the systematic uncertainty. In appendix A the con¬ 
tributions from each source of systematic uncertainty, described in section 7, to the total 
systematic uncertainty in each bin of the A|^| and A|y| distributions are provided. The 
statistical correlations between the different bins of the distributions are also given. 

The results for the inclusive lepton-based asymmetry Aq and the tt asymmetry Aq 
after corrections for detector and resolution effects are shown in table 6. The values in 
the ee, e/r and py, channels as well as for their combination are presented, together with 
statistical and systematic uncertainties. 

Detailed information about the combination of the inclusive values is given in table 7. 
The combination probabilities are 21% and 81% for and Aq respectively, demonstrating 
the compatibility of the measurements in the three channels (ee, ep and /ip). The weight 
of each channel in the combination is also reported in table 7. The ep channel dominates 
the combination, reflecting the larger data statistics compared to that of the ee and pp 
channels. 

The inclusive measurements after the detector and resolution effects corrections can be 
compared with the state-of-the-art theoretical predictions calculated at NLO QCD, includ¬ 
ing the electromagnetic and weak-interaction corrections [10]: Aq = 0.0070 ±0.0003 (scale) 
and Aq = 0.0123 ± 0.0005 (scale). In figure 7 the measured values of Aq and Aq are 
compared to these predictions and PoWHEG-hvq ± PYTHIA6 predictions. In the figure, el- 
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Bin of A\r]\ 


1 dc7 


(7 dA|r;| 


(± stat. ± syst.) 


[-3.00,-2.00] 

0.0440 

± 

0.0077 

± 

0.0025 

[-2.00,-1.67] 

0.106 

± 

0.011 

± 

0.004 

[-1.67,-1.33] 

0.126 

± 

0.011 

± 

0.005 

[-1.33,-1.00] 

0.164 

± 

0.012 

± 

0.004 

[-1.00,-0.67] 

0.245 

± 

0.013 

± 

0.004 

[-0.67,-0.33] 

0.314 

± 

0.015 

± 

0.007 

[-0.33,0.00] 

0.400 

± 

0.016 

± 

0.004 

[0.00,0.33] 

0.392 

± 

0.016 

± 

0.004 

[0.33,0.67] 

0.349 

± 

0.015 

± 

0.009 

[0.67,1.00] 

0.244 

± 

0.013 

± 

0.010 

[1.00,1.33] 

0.209 

± 

0.013 

± 

0.005 

[1.33,1.67] 

0.129 

± 

0.011 

± 

0.003 

[1.67,2.00] 

0.0815 

± 

0.0093 

± 

0.0028 

[2.00, 3.00] 

0.0361 

± 

0.0076 

± 

0.0035 


Table 4: Normalized differential cross-sections for A|r/| in the e/r channel presented with 
statistical and systematic uncertainties. 


Bin of A y 

ld%y\ (± stat - ± sy st -) 

[-5.00,-0.71] 

[-0.71,0.00] 

[0.00,0.71] 

[0.71,5.00] 

0.0435 ± 0.0029 ± 0.0017 

0.437 ± 0.016 ± 0.013 

0.420 ± 0.025 ± 0.015 

0.0470 ± 0.0032 ± 0.0024 


Table 5: Normalized differential cross-sections for A|y| in the ey channel presented with 
statistical and systematic uncertainties. 


Channel 

4? 

4 

ee 

0.101 ± 0.052 ± 0.021 

0.025 ± 0.069 ± 0.027 

ey 

0.009 ± 0.019 ± 0.009 

0.007 ± 0.032 ± 0.018 

MM 

0.047 ± 0.030 ± 0.012 

0.043 ± 0.045 ± 0.013 

Combined 

0.024 ±0.015 ±0.009 

0.021 ±0.025 ±0.017 

SM, NLO QCD+EW [10] 

0.0070 ± 0.0003 (scale) 

0.0123 ± 0.0005 (scale) 


Table 6: Results for the lepton-based asymmetry Aq and the tt asymmetry Aq after 
correcting for detector, resolution and acceptance effects. The values in the ee, e/j, and /ry 
channels as well as the combined value are presented with their statistical and systematic 
uncertainties. 

lipses corresponding to lcr and 2a combined statistical and systematic uncertainties of the 
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aU 

A C 

Mg 

x 2 

Probability (in %) 

Weights (ee/efj,/fifj, in %) 

3.1 

21 

7/68/25 

0.4 

81 

9/57/34 


Table 7: Information about the combination of the three channels using the best linear 
unbiased estimator method: x 2 and probability of the combination, as well as the weight 
of each channel. 
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Figure 7: Comparison of the inclusive Aq and Aq measurement values to the theory 
predictions (SM NLO QCD+EW prediction [10] and the prediction of the PoWHEG-hvq 
+ PYTHIA6 generator). Ellipses corresponding to la and 2a combined statistical and 
systematic uncertainties of the measurement, including the correlation between Aq and 
Aq, are also shown. 
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measurement, including the correlation between Aq and A g, are also shown. The statis¬ 
tical correlation between Aq and A g is evaluated to be 37+5% using pseudo-experiments 
based on simulation. The systematic uncertainties are treated as 100% correlated. The 
resulting correlation between Aq and A g is about 55%. The measured values are both 
consistent with the theory predictions within the uncertainties. The measured A g values 
are consistent with but less precise than measurements in the single-lepton decay channel 
by the ATLAS [18] and CMS [19] collaborations. The measurements of Aq and Mg are also 
consistent with the CMS collaboration measurements in the dilepton decay channel [20]. 

The inclusive measurement of Aq and Aq is furthermore compared to two models 
of physics beyond the Standard Model (BSM) [9] that could be invoked to explain an 
anomalous forward-backward asymmetry at the Tevatron, such as reported by the CDF 
experiment [24], Two models with a new colour octet particle exchanged in the s-channel 
are considered. In the model with the light octet, the new particle mass is below the tt 
production threshold. The model with the heavy octet uses the octet mass beyond the reach 
of the LHC. The new particles would not be visible as resonances in the m t i spectrum at the 
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Figure 8: Comparison of the measured inclusive Aq and Aq values to two benchmark 
BSM models, one a light octet with mass below tt production threshold (left) and one with 
a heavy octet with mass beyond LHC reach (right), for various couplings as described in 
the legend. 


Tevatron or at the LHC. The light octet is assumed to have a mass of m = 250 GeV and a 
width of r = 0.2m. For the heavy octet, the corrections to tt production are independent of 
the mass but instead depend on the ratio of coupling to mass, which is assumed to be 1/TeV. 
In figure 8 the measured Aq and Aq values are compared to the light (figure 8(a)) and 
heavy (figure 8(b)) colour octet model predictions in order to assess whether any of the BSM 
predictions can be excluded. Models with left-handed, right-handed and axial coupling to 
the up, down and top quarks are shown. The considered couplings to the quarks are such 
that the global fit to tt observables at the Tevatron and the LHC, including total cross- 
sections, various asymmetries, the top polarisation and spin correlations, is consistent with 
the measurements within two standard deviations [9]. The LHC asymmetry measurements 
in the dilepton decay channel are excluded from the fit. While the models span a sizeable 
range of values in the Aq and Aq plane in figure 8, their predictions are consistent with the 
measured value within the present uncertainties. Thus the potential BSM contributions 
cannot be excluded beyond the reach of the previous Tevatron and LHC measurements. 
Future Aq and Aq measurements with a larger dataset could however further constrain the 
allowed couplings of the colour octet models if both statistical and systematic uncertainties 
can be reduced further. 

9 Conclusion 

Measurements of the tt charge asymmetry in the dilepton channel are presented. The mea¬ 
surements are performed using data corresponding to an integrated luminosity of 4.6 fb -1 
of pp collisions at y/s = 7 TeV collected by the ATLAS detector at the LHC. Selected 
events are required to have exactly two charged leptons (electron or muon), large miss- 
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ing transverse momentum and at least two jets. Both the lepton-based asymmetry Aq 
and the tt asymmetry Aq are extracted in three channels: ee, e/x and /x/x. The mea¬ 
surement of Aq requires the kinematic reconstruction of the tt system, which is per¬ 
formed using the neutrino weighting technique. Agreement between predictions and data 
is checked after selection and kinematic reconstruction. Good agreement is obtained for 
all the kinematic observables studied. The ££ A\tj\ and tt A\y\ distributions and inclu¬ 
sive asymmetries are corrected for detector and acceptance effects. Corrections are ap¬ 
plied using bin-by-bin corrections for Aq and fully bayesian unfolding for Aq. The dis¬ 
tributions of lepton A\r]\ and ti A\y\ after the detector smearing corrections are provided 
for the e/x channel. Good agreement between the corrected values and predictions of the 
Monte Carlo generator models is observed in these distributions. The combined values of 
lepton-based inclusive asymmetry 4* and tt inclusive asymmetry Aq are measured to be 
A% = 0.024T0.015 (stat.)±0.009 (syst.) and Ag = 0.021±0.025 (stat.)±0.017 (syst.). The 
measured values are in agreement with previous LHC measurements and with the Standard 
Model prediction [10]: Aq = 0.0070±0.0003 (scale) and A(? = 0.0123±0.0005 (scale). The 
measurements are limited by statistical uncertainties. The predictions of benchmark light 
and heavy colour octet models with parameters selected such that the models are consistent 
with previous LHC and Tevatron data [9] are found to be consistent with the measured 
asymmetries. 
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A Additional tables 


Additional information about the normalized differential cross-sections in the e/x channel 
are provided in this appendix. 

The detail of the systematic uncertainties in each bin of the distributions are reported 
in tables 8 and 9. 


Bin of A|? 7 | 

[-3.,-2.] 

[-2.. -1.67] 

[-1.67, -1.33] 

[-1.33, -1.] 

[-1., -0.67] 

[-0.67, -0.33] 

[-0.33,0.] 

Central value 

Stat. 

0.0440 

0.0077 

0.106 

0.011 

0.126 

0.011 

0.164 

0.012 

0.245 

0.013 

0.314 

0.015 

0.400 

0.016 

Lepton reconstruction 

0.0023 

0.003 

0.005 

0.003 

0.003 

0.005 

0.003 

Jet reconstruction 

0.0005 

0.001 

- 

0.002 

- 

0.002 

0.001 

r-imiss 

il/rji 

- 

- 

- 

- 

- 

- 

- 

Signal modelling 

n/r 

n/r 

n/r 

n/r 

n/r 

n/r 

n/r 

PDF 

n/r 

n/r 

n/r 

n/r 

n/r 

n/r 

n/r 

Background 

- 

0.001 

0.001 

- 

0.001 

0.002 

0.002 

NP and fake 

0.0005 

0.002 

0.001 

- 

0.001 

0.004 

0.001 

Total systematics 

Total uncertainty 

0.0025 

0.0081 

0.004 

0.011 

0.005 

0.012 

0.004 

0.013 

0.004 

0.014 

0.007 

0.016 

0.004 

0.017 


Bin of A|? 7 | 

[0., 0.33] 

[0.33,0.67] 

[0.67,1.] 

[1., 1.33] 

[1.33,1.67] 

[1.67,2.] 

[2., 3.] 

Central value 

0.392 

0.349 

0.244 

0.209 

0.129 

0.0815 

0.0361 

Stat. 

0.016 

0.015 

0.013 

0.013 

0.011 

0.0093 

0.0076 

Lepton reconstruction 

0.003 

0.009 

0.009 

0.003 

0.003 

0.0019 

0.0032 

Jet reconstruction 

0.001 

0.001 

0.002 

0.002 

- 

0.0012 

0.0003 

T-tmiss 

- 

- 

- 

- 

- 

- 

- 

Signal modelling 

n/r 

n/r 

n/r 

n/r 

n/r 

n/r 

n/r 

PDF 

n/r 

n/r 

n/r 

n/r 

n/r 

n/r 

n/r 

Background 

0.002 

- 

- 

0.001 

0.001 

0.0009 

- 

NP and fake 

0.001 

0.001 

0.002 

0.003 

0.001 

0.0012 

0.0013 

Total systematics 

0.004 

0.009 

0.010 

0.005 

0.003 

0.0028 

0.0035 

Total uncertainty 

0.017 

0.018 

0.017 

0.014 

0.011 

0.0097 

0.0084 


Table 8: Systematic uncertainties in each bin of the A|? 7 | distribution in the e/x channel. 
Hyphens are used when the uncertainties are lower than 0.0005. The signal modelling 
and the PDF uncertainty, (labeled as n/r) are limited by the statistical fluctuations in the 
simulated samples and are thus not reported in the table. 

The statistical correlations between the different bins of each distribution are reported 
in tables 10 and 11. They were estimated using bootstrapping. 
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Bin of A|j/| 

[-5.00,-0.71] 

[-0.71,0.00] 

[0.00,0.71] 

[0.71,5.00] 

Central value 

0.0435 

0.437 

0.420 

0.0470 

Stat. 

0.0029 

0.016 

0.025 

0.0032 

Lepton reconstruction 

0.0009 

0.006 

0.010 

0.0016 

Jet reconstruction 

0.0007 

0.011 

0.009 

0.0012 

rpmiss 

Jl /'p 

0.0002 

0.003 

0.007 

0.0009 

Signal modelling 

n/r 

n/r 

n/r 

n/r 

PDF 

n/r 

n/r 

n/r 

n/r 

Background 

0.0003 

0.001 

0.001 

0.0002 

NP and fake 

0.0011 

0.001 

0.004 

0.0007 

Total systematics 

0.0017 

0.013 

0.015 

0.0024 

Total uncertainty 

0.0034 

0.021 

0.029 

0.0040 


Table 9: Systematic uncertainties in each bin of the A|y| distribution in the e/x channel. 
The signal modelling and the PDF uncertainty, (labeled as n/r) are limited by the statistical 
fluctuations in the simulated samples and are thus not reported in the table. 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

1 

+ 1.00 

-0.49 

-0.02 

-0.03 

-0.02 

+0.01 

-0.03 

+0.02 

+0.06 

-0.06 

-0.01 

+0.03 

-0.02 

+0.03 

2 

-0.49 

+ 1.00 

-0.50 

-0.04 

+0.02 

+0.00 

+0.00 

+0.01 

-0.04 

+0.05 

-0.02 

-0.01 

+0.01 

-0.03 

3 

-0.02 

-0.50 

+ 1.00 

-0.47 

+0.00 

+0.02 

-0.03 

-0.02 

-0.03 

+0.03 

+0.03 

+0.01 

-0.02 

-0.00 

4 

-0.03 

-0.04 

-0.47 

+1.00 

-0.49 

-0.04 

-0.00 

+0.04 

+0.03 

-0.05 

+0.03 

-0.03 

+0.03 

-0.01 

5 

-0.02 

+0.02 

+0.00 

-0.49 

+ 1.00 

-0.52 

+0.06 

-0.04 

-0.01 

+0.03 

-0.05 

+0.00 

-0.04 

+0.06 

6 

+0.01 

+0.00 

+0.02 

-0.04 

-0.52 

+ 1.00 

-0.54 

+0.04 

-0.02 

+0.01 

+0.01 

-0.00 

+0.02 

-0.01 

7 

-0.03 

+0.00 

-0.03 

-0.00 

+0.06 

-0.54 

+ 1.00 

-0.54 

-0.00 

-0.00 

-0.00 

+0.04 

-0.02 

+0.01 

8 

+0.02 

+0.01 

-0.02 

+0.04 

-0.04 

+0.04 

-0.54 

+1.00 

-0.53 

+0.02 

+0.02 

-0.05 

+0.00 

-0.02 

9 

+0.06 

-0.04 

-0.03 

+0.03 

-0.01 

-0.02 

-0.00 

-0.53 

+1.00 

-0.52 

+0.01 

+0.00 

+0.03 

-0.02 

10 

-0.06 

+0.05 

+0.03 

-0.05 

+0.03 

+0.01 

-0.00 

+0.02 

-0.52 

+ 1.00 

-0.52 

+0.02 

-0.02 

+0.01 

11 

-0.01 

-0.02 

+0.03 

+0.03 

-0.05 

+0.01 

-0.00 

+0.02 

+0.01 

-0.52 

+ 1.00 

-0.49 

+0.03 

-0.03 

12 

+0.03 

-0.01 

+0.01 

-0.03 

+0.00 

-0.00 

+0.04 

-0.05 

+0.00 

+0.02 

-0.49 

+ 1.00 

-0.56 

+0.02 

13 

-0.02 

+0.01 

-0.02 

+0.03 

-0.04 

+0.02 

-0.02 

+0.00 

+0.03 

-0.02 

+0.03 

-0.56 

+ 1.00 

-0.54 

14 

+0.03 

-0.03 

-0.00 

-0.01 

+0.06 

-0.01 

+0.01 

-0.02 

-0.02 

+0.01 

-0.03 

+0.02 

-0.54 

+ 1.00 


Table 10: Statistical bin-bin correlations within the A\r/\ distribution in the ejx channel. 
The bin numbers are used instead of the bin boundaries. Bin 1 corresponds to [—3.00, —2.00] 
and bin 14 to [2.00,3.00]. 
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1 

2 

3 

4 

1 

+1.00 

-0.63 

+0.38 

-0.38 

2 

-0.63 

+1.00 

-0.79 

+0.37 

3 

+0.38 

-0.79 

+1.00 

-0.61 

4 

-0.38 

+0.37 

-0.61 

+1.00 


Table 11: Statistical bin-bin correlations within the A|y| distribution in the e/i channel. 
The bin numbers are used instead of the bin boundaries. Bin 1 corresponds to [—5.00, —0.71] 
and bin 4 to [0.71, 5.00]. 
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